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Abstract 

We study n-point boundary correlation functions in Timelike Boundary Li- 
ouville theory, relevant for open string multiproduction by a decaying unstable 
D-brane. We give an exact result for the one-point function of the tachyon ver- 
tex operator and show that it is consistent with a previously proposed relation 
to a conserved charge in string theory. We also discuss when the one-point am- 
plitude vanishes. Using a straightforward perturbative expansion, we find an 
explicit expression for a tachyon n-point amplitude for all n, however the result 
is still a toy model. The calculation uses a new asymptotic approximation for 
Toeplitz determinants, derived by relating the system to a Dyson gas at finite 
temperature. 
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1 Introduction and Summary 



String theory contains D-branes of opposite charge, so one should be able to under- 
stand their annihilation process. A related problem is the decay of a single unstable 
brane, such as a D-brane in bosonic string theory A simple model for the D-brane 
decay describes a process starting from the infinite past, involving a spatially homo- 
geneous tachyon field rolling towards the true minimum of its effective potential [1,2]. 
A basic open problem is to calculate tree-level string scattering amplitudes in the 
rolling tachyon background, corresponding to production of multiple closed or open 
strings by the decaying brane. There are both conceptual and technical aspects to 
this problem. Because the background is time dependent, there are different ways to 
define the notion of vacuum and asymptotic states. A technical framework for the 
bosonic homogeneous brane decay is the timelike boundary Liouville theory (TBL) 
coupled to 25 free massless spacelike bosons, and the problem of computing n-point 
correlation functions [3]. Calculations are difficult since they involve complicated 
coupled integrals and/or nonintuitive analytic continuations. 

In this paper we focus on calculating boundary n-point functions in TBL. The 
two-point function, associated to the rate of open string pair production by a decay- 
ing brane, has been investigated before [3,4], and also in a curved spacetime (AdSs) 
in [5]. (Other string production work is found in [6-16].) A simple toy model is 
obtained by moving to the minisuperspace approximation, where strings are point- 
like, and the problem reduces to a relatively simple quantum mechanical scattering 
problem. Returning back to the original setup, the standard prescription is to start 
from spacelike boundary Liouville theory (SBL), where the two-point and three-point 
functions have known well-defined analytic expressions [17-19], and then continue 
to the timelike theory by b — > i, <fi — > iX°. However, the continuation must involve 
a prescription to avoid the accumulation of an infinite number of poles and zeroes 
which would render the answer ill defined. One way to motivate a prescription is 
by aiming to make contact with the minisuperspace analysis. This procedure gives 
a physically pleasant answer, exponentially suppressed open string pair production 
at high energies. However, some doubt remains, as the prescription for the analytic 
continuation was not unique and some of the steps involved are rather indirect. It 
is desirable to pursue alternative approaches, they may give further support to the 
previous analysis or lead to other reasonable prescriptions. Moreover, the previous 
method is difficult to extend beyond the two-point function. 

An alternative method to compute correlation functions in TBL was given in 
[15]. Instead of indirect arguments, the method [15] is based on a straightforward 
perturbative expansion, and the observation that Random Matrix Theory (RMT) [20] 
techniques become applicable to the ensuing integrals. This method was successfully 
applied to compute the bulk-boundary function [15, 16]. On the other hand, the 
same problem was also considered by Liouville theory methods. The bulk-boundary 
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function was calculated in spacelike Liouville theory in [21]. Ref. [22] then investigated 
the analytic continuation from spacelike Liouville theory to timelike theory, and found 
a result for the bulk-boundary function which is similar to that of [15]. 

In the method of [15], correlation functions are related to expectation values of 
periodic functions (Fisher-Hartwig symbols) in the circular ensemble of unitary ma- 
trices (CUE), also equivalent to Toeplitz determinants of Fourier coefficients. This 
observation was extended to n-point functions and superstrings in [16]. Alternatively, 
the rz-point functions can be related to thermal expectation values in a classical log 
gas of unit charges in two dimensions, e.g., the Dyson gas. In [15], this observation 
was made at a formal level, while the problem of actually finding explicit answers for 
the correlation functions still remained. 

In the present paper, we use the interpretation of the correlations functions as 
thermal Dyson gas expectation values, and then use physical insight to find analytic 
expressions. We are able to derive an expression for an n-point amplitude. The virtue 
of our approach is that it is relatively straightforward, and it is powerful enough 
to for a first time yield an analytic expression for an n-point amplitude for all n. 
The downside is that at the moment we do not have quantitative control of our 
approximations by the time we compute the amplitude. Consequently, we do not 
yet know how to compare the result with the previous one for the open string pair 
creation amplitude. Nevertheless, we consider the techniques that we have developed 
to be a step forward towards full control of the scattering problem. 

The paper is organized as follows. Section [2] begins with a review of some facts 
of TBL. We present some preliminary calculations and discuss the one-point function 
and the vanishing one-point amplitude. Next, we present a contour integration trick 
which is powerful in summing the series expansion for the correlation functions. We 
then use the one-point function formula to test a recently proposed master formula [23] 
for a string theoretic definition of a conserved charge. Section [3] reviews the relation 
of TBL to Dyson gas at finite temperature. In Section HI we use this connection to 
derive an approximation for the integrals which appear as coefficients in the series 
expansion of a n-point amplitude, then use the approximation for the coefficients 
and the contour integration trick of Section [2] to derive a toy model result for the 
amplitudes. Some calculational details are left in Appendices |Aj |Bl and O 

2 Boundary amplitudes in timelike Liouville the- 
ory 

Let us first review some facts to identify the problem of interest. Full scattering 
amplitudes in bosonic string theory involve contributions from the timelike X° and 
the 25 spacelike directions X = (X 1 ). However, as discussed in [7], one can simplify 
the calculations by adopting a gauge where the string vertex operators factorize into 
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a form 

V = e^ x ° V sp (X) , (1) 

so that all dependence on X° is in the simple exponential factor, while V sp contains 
the more complicated polarization tensor factor and only depends on the spacelike 
directions X. For a homogeneous rolling tachyon background depending only on 
X°, all the complications arise from contractions in the X° direction between the 
background and the vertex operators, while contractions in the spatial directions 
give a simple contribution. Correspondingly, the n-point correlation functions in 
the homogeneous rolling tachyon background factorize into a product of an n-point 
function of e tuJ aX°(T a ) ( w ]j ere label a — 1, ... ,n) in TBL, and an n-point function 
of V sp (k a ; X{r a )) in the theory of free spacelike bosons, 

In \ In 

JP^A m Vsp (k a ;X(r a )) 

\a=l I TBL \o=l 

= e -iE a ^m e ^o (ra) \ Ffree [(4) ; ( Ta) ] ; (2 ) 

\a=l / TBL 

where we separated the spacelike zero modes. The on-shell conditions k 2 a = —uj^+kl = 
—m 2 a can be satisfied for a range of values of u a . The problem of interest is to calculate 
n-point functions in TBL for generic u a . We will also try to compute the full scattering 
amplitude for n open string tachyons. 
The action of TBL is 

S TBL = -— [ dX°dX° + A / dte x ° . (3) 




2tt 



disk 



Eventually we will be interested in the open string n-point tachyon amplitude 

A n ^i,ki,...;oo n ,k n ) = [ (Pxe-^^ I f[pLF {ree [(k a );(r a )} (4) 

J J a=l 71 



/IL 
VX°e~ STBL Y[ e iuJaX ° iTa) 
a=l 



at tree level, where the momenta k a are in the spacelike directions of the decay- 
ing jo- dimensional brane, r a denote points on the boundary of the disk (unit cir- 
cle) a For tachyons the contribution from the spacelike directions (with divergent 
self-contractions removed) is 

F bee {(k a ); (r a )} = J] \e iT ' - ^f^» , (5) 

a<b 



1 We could use the conformal Killing group (CKG) PSL(2, R) to fix three of the vertex operator 
coordinates t , but we have chosen to leave them unfixed and average over the locations. 
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with the on-shell condition k\ = —ul+kl = 1. The conservation of spatial momentum 
has been discussed, e.g., in [15]. As discussed in the Introduction, different approaches 
have been used for the calculation. We will follow the approach of [1, 15] and first 
expand A n as a power series, in powers of the boundary interaction. We also separate 
out the overall zero mode a; dependence, so A n becomes 



(2=1 

x/n| 



dx°e ix °^=^ a 

N=0 



(-2ir\e x °) 



x ,0 (t N ) 



n 



a=l 



,iLO a X'°{T a ) 



(6) 



After the Wick contractions and substituting the Green's functions, the amplitude 
takes the form of a power series of coupled integrals. 
The amplitude A n then becomes 



dx° exp 



a=l 



A n (2n\e x ) , where 



N=0 

where we have adopted the notation 

z = 2n\e x ° 

and the integrals 



(7) 



(8) 



X 



2tt ax 2tt 

i=l a=l 



n 

l<i<j<N 



ffii I ^ 



nn 

i=l a=l 



l<a<fe<n 



, (9) 



which include the spacelike contribution F. In order to do the sum over N we need 
to work out the tj integrals for arbitrary N. When calculating the integrals it is often 
useful to assume that £ a are positive real numbers and continue to imaginary £ a , i.e., 
to real energies cu a , in the end. This is not problematic since the ti integrals converge 
for Re£ a > —1/2 and thus define an analytic function of £ a in this region. 



2.1 Some preliminary considerations 

The simplest case to consider is n — 1, the one-point boundary amplitude. Invariance 
under translation requires the one-point function to vanish unless the operator at the 
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boundary has zero conformal weight, rendering the case trivial. However, it turns out 
that some calculations will be useful for the nontrivial case n > 1. It is also known 
that in a noncompact conformal field theory (CFT) integrated one-point functions 
can be nonzero [23,24]. Ref. [23] proposed a relation between a one-point function 
and a spacetime boundary term. In our case, we can use a TBL one-point function 
as a check of the master formula in [23], and find it to be consistent. 

Let us postpone other discussions for a moment and just focus on a straightforward 
calculation. We consider the series that appears in ©, with n = 1, 



A 1 (z) = J2(-z) N -h(N), 



(10) 



N=0 



where now 



L(N) 



1 

M 
l 

m 



2tt 



dr 
2^ 



N 



i=l 



dtj 
2?r' 



II 



iti I 2i(j 



n 

l<i<j<N 



g^i I ^ 



i=l 



2tt 



n 



iti 1 2 



X<i<j<N 



(12) 



Here we denoted £ = iaj, where u is the energy of the open string^] It is interesting 
to note that the integrand is independent of r, the coordinate of the vertex operator 
at the boundary, so that the r integral is trivial. In other words, the integrand is 
invariant under translations along the boundary, independently of £ — iu. However, 
the total one-point function also contains the contribution from the spacelike direc- 
tions with a 8o,k factor, which along with the on-shell condition will constrain u. But 
let us focus back to the properties of the series ffTUj) . 

The same series has been considered in the context of a general bulk-boundary 
amplitude, which has been calculated in closed form in [15,16]. The bulk-boundary 
amplitude involves 



/oo 
dx o e ^o+^°J2(-z) N I luJo (N) 
N=0 



(13) 



where u c is the energy of the bulk operator exp{iu> c X°(z, z)} and I iuJo {N) is the 
integral (fTTj) evaluated at £ = iu , where u Q is the energy of the boundary operator Jf| 
First, the integral evaluates to the relatively simple expression 



N 



im = n 



r(j)T(j + 20 G(f + l) 2 G(N + 2£ + l)G(N + 1) 



r( J + o 2 G(2£ + i; 



G(N + z + iy 



(14) 



2 Note that (after removing the self-contractions in the spacelike directions) Ff roc = 1. 
3 The one-point amplitude is formally the limit lo c — > of the bulk-boundary amplitude (|13| . We 
also omitted a ^-function term [see Subsection 12.2] . 
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where G is the Barnes G function. After converting to integral representation of the 
T functions, the sum over N in (imp can be done [15, 16], leading to the result 



Ai +1 (uj c ,u ) 



-Z7T- 



(2ttA) 



-i(w c +LU ) 



sinh Ti (u c + lu 



exp 



dt(l - e 



-iu) t\2 



J 2t(l-cosht) 



3 i(aj c +u; )i\ 



• (15) 



We would first like to point out an interesting feature, which was not investigated 
in [15, 16]. Let us write it in terms of the Barnes G functions, using the integral 
representation 

-zt 



iogG(* + i; 

We find 
where 



dt 
— < 

t 



zz-r 



l-e 



-t 



+ 



1 



:i-e- 



Ai+i(uj c , uj ) = -in , ■ -J iuJo (i(u c + uj )) 



sinh 7t(uj c + u 



Re(z) > -1 . 

(16) 

(17) 



Jt(s) 



G(e + l) 2 G{2j - s + l)G(-s + 1) 



G(2f+1) G(e-s + l) 2 
The asymptotic behavior [15,16] follows easily from f|T8l) . 



(19) 



by using the asymptotic series of the Barnes G function 

logG(z + l)=z 2 Qlogz-|) +^lo g 2vr-llog^ + C , (-l) + 0(lA 2 ) • (20) 

An interesting feature is that ([TBI) is a natural continuation of ()14p to noninteger 
values, replacing iV — > —s, but (TH1) was the A^th coefficient in the series ( 1131) . while 
([TH[) is essentially the sumQ We will show how coefficients convert to the sum in 
the next Subsection 12. 2[ by a new contour integral trick which also allows a more 
controlled investigation of the convergence of the series (fl3[) . The other benefit of 
the calculation is that it can also be applied to n-point amplitudes. But let us first 
continue with the one-point function. 

As seen by comparing f[T0[) and f[T3l) . we can formally use the result (|17|) to obtain 
a formula for the Fourier transform of ([TDl) by setting uj c — and io — u, giving 



Ai(u) 



A 1+1 (0,u) 
(2ttA)- 



-ITX- 



sinh tiuj 



exp 



)(1 — cosut) 



00 , (1 - e-*"* 

dt 7 — 

o t(l — cosht) 



(27rA)-^r(icj 



G(zcj + 1) 3 G(2-z^) 
' G(2zo; + l) 



(21) 



4 A similar observation has been made in the case of bulk amplitudes in spacelike Liouville theory 
in [25]. 
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Notice that we will carefully rederive this formula in the next subsection. The sin- 
gularities and zeroes of this function are listed in Appendix In particular, the 
zeroes are located at the imaginary axis, at u = in, where n is an integer, except at 
n = 0, ±1. Consider then the full one-point tachyon amplitude (the n — 1 case of 

Ai(u) = 6 £ J dx° exp(iLux°)Ai(2ir\e x °) = S 0)k Ai{uj) . (22) 

The momentum conservation condition k = along with u 2 = — 1 + k 2 demands 
uj = ±z so that the amplitude involves the operator exp(=|=X ). The result is 



AM=5 ^( 7 rA) ±1 . 123) 



Note that the choice u = — i is related to the disk partition function by 

M<» = -i) = ~ / ^l; W* ) = ^ , («) 

where ^disk,A(a ;0 ) = Aq(x°) = 1/(1 + 2nXe x °). Conversely, for uj ^ ±i, the on-shell 
condition requires k ^ so that the one-point amplitude vanishes. Even though the 
amplitude vanishes for generic u, the expression (12T|) will be met again in the context 
of higher point amplitudes. It will be interesting to know its asymptotic behavior in 
the limit |o;| — > oo. It can be calculated to arbitrary order by using the asymptotics 
of Barnes G (1201) . The leading terms are 



-IUJ 



A I S ■ ( 27rA )" 

Auoj) = —m exp 

sinh huj 



uj 2 ( ySgn(Reu;) + 21og2j - - log (iu) 



-sgn (Reu;) + -log2 + 3C / (-i; 



[1 + 0{uj~ 2 )} , (25) 



where arg(u;) 7^ ±7r/2. 



2.2 A contour integral method 

Next we calculate the integrated amplitude using a contour integration trick which 
allows us to sum the series over N in ([7j) and analytically continue the resulting 
amplitude to the region where the defining sum is not convergent. The essential 
required feature of the coefficients I%(N) is that they should not diverge too fast for 
large N. For concreteness and simplicity we will first consider the series ffTOl) and f[T5|) . 
However, our method can also be applied to higher point functions as we will discuss 
in Section HI More precisely, the calculation can be generalized to the case of the n- 
point amplitude © if we use a suitable approximate form for the integral coefficients 
l£ 1 ,...,£ n (N) . As the contour integration method enables us to control the convergence 
of the sum and the integral it is more rigorous than the original calculation in [15,16]. 



8 



We begin by studying the analytic structure of J^(s) of (1TB]) and the asymptotics of 
I%{N) for large N. We will first consider the case where £ is real and positive. Recall 
the continuation of the coefficient formula (j!4p to noninteger values of N = — s, given 
by ([IS]) . From the asymptotic formula of the Barnes G funtion (120]) it immediately 
follows that J^(s) has a powerlike behavior for large s, 



G(e + i) 2 2 

fl j G(2£ + l) 1 j 



s 



args^O. (26) 



In addition, since G(z + 1) is an entire function with zeroes at z — —1, —2, . . ., the 
poles of Jf(s) are located a 

Thus in the region \z\ < 1, where the sum in ( [13]) converges, the asymptotic 
behavior of in f[2"B"]) enables us to write the sum as 

AT=0 ^ C 

where the contour C wraps around the negative real s axis as depicted in Figure [TJ 
picking up the residues at the poles of 1/ sin(7rs) at s — 0, —1, —2, . . . which produce 
the terms in the series. Note that the zeroes of G(l — s) in J^{s) cancel the poles 
of l/sin(7rs) for s = 1,2,3,.... Since l/sin(7rs) vanishes exponentially for large 
imaginary s, we may deform the contour (keeping \z\ < 1) in (|27|) to 

A 1 (z) = — 1?—J 6 (s)d8, (28) 

Jj-ioo SUMS 

where < 7 < £ + 1. This integral converges everywhere except for negative real z 
(if the principal branch of z~ s with | argz| < 7r is used) and thus defines the analytic 
continuation of A\{z) to \z\ > 1, | argz| < n. Moreover, for \z\ > 1 we can continue 
to deform the contour to 

A 1 {z) = ±-l ^Js(s)ds, (29) 
2m J c , sm7rs 

where C wraps around the positive real s axis as shown in Figure [TJ The integral 
is convergent for all \z\ > 1 so there are no singularities in this region but a loga- 
rithmic branch cut ending at z = 00 which arises from the factor z~ s . The residue 
contributions at the poles of J^(s) at s = £ + 1, £ + 2, . . . give the \ jz expansion 



Ai(z) = (Q + ATogz)^ 1 



l + O 



(30) 



where the constants C^D^ can be calculated using (TT8]) . 



3 For £ = 1,2,... the poles are found at s = 2£, 2£ + 1, . 
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Figure 1: The different integration contours on the s plane that define the analytic 
continuation of A\(z) for all values of z. Integration over the contours C, 7 -MR, and 
C, converge for \z\ < 1, |arg2;| < ir, and \z\ > 1, respectively. The x's and the o's 
denote the poles of Jg(s)/ sin7rs. 

To summarize, from the different integral representations (j27j) . (|28|) . and (1291) it 
follows that the only singular points of Ai(z) are z — — 1 and z = 00. In particular, 
on the integration path in the one-point amplitude 

/oo 
dx^^A^x ) (31) 
-00 

i.e., z = 2n\e x ° = ... 00, A\(z) has no singularities. Using the series in (T27|) and in 
(pOj) we see that the integrand vanishes exponentially 



x u — »oo 

^ ~ (32) 



for large ±a;° so the integral over x in fl3Tl) is convergent. Moreover, note that 
inserting the definition of z in (|28[) we find 



A(x u )=/ n{2 * Xr ' e ~ Ms)ds (33) 
which defines the inverse of the (bilateral) Laplace transform. The inverse relation 
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then gives the master formula for the one-point amplitude in terms of J^{s), 



Ms) = I dx°e sxU Ai(x°) = ^ 27rX) ' J*(s) . (34) 

The steps from ( 1271) to (134")) show how the analytic continuation of the coefficients 
of the series ends up as its sum. From the asymptotics of Ai(x°) we see that (13~4"1) 
converges for < s < £ + 1, as expected from the positions of poles of J{(s)/ sin7rs 
[and the choice of 7 in (1281)]. In particular, 

A(0 = As = = = (2ttA) e r(e) G(2e + 1) (35) 



reproducing the result ( 1211) above. 

In the end, we want to continue the result (1331) for the integrated amplitude for 
imaginary £ = iuo. For imaginary £ the above analysis is not essentially changed: the 
poles of J%(s) move to s = iu + 1, iuj + 2, . . ., but still lie to the right of the imaginary 
axis, so that A~i(x°) vanishes exponentially Ai(x°) ~ e~ x ° for a; — > 00. However, 
after inserting s = iu in ( 13~4"1) the convergence in the opposite direction x° — * — 00 is 
lost. We find instead 

e^AAx ) . ~ e^ x '° (36) 
which signals the presence of a 5 function. Indeed, the integral can be interpreted afl 

G(l+iu) 3 G(2-iu) 



G(2iu + 1) 



- wrw-* )) C(1+ G g^ 2 1) - M . m 

where the changes the value of slightly to that direction where the x° integral is 
convergent. 



2.3 The one-point function as a boundary term in spacetime 

As discussed in [23], one difference between CFTs in compact and noncompact target 
spacetimes is that in the latter case boundary terms can spoil the holomorphicity of 
the stress tensor. This modifies its OPE with other operators, and lead [23] to derive 
a master formula relating the one-point function (on a sphere or at the boundary of 
a disk) to a boundary term in spacetime, so as to give a string theoretic definition for 
a conserved charge, as an extension from field theory. For a disk one-point function, 

6 The result can be checked explicitly by writing Ai(x°) = AiX^ )^^ + [-4i(:e ) — -^i( :c0 )| w=0 ] 
where the first term is simple to integrate and the latter does not contribute to the singularity. 
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the master formula is 



(0(z,z)) 



M I d D x d, 



d 2 z'e 2 ^ z '~ z,) 



D 2 



+ 



z' + z 
2z 

z' + z 
2z 



(z' - z)(dX^(z',z r )0(z,z)) 



D 2 



(z'-z){dX»(z>,z')0(z,z)) D2 



(38) 



where J\f is a normalization factor, the metric on the disk is ds 2 = e 2u) ( z ' z ^dzdz' and 
0(z, z) is a local boundary operator in the CFT with D- dimensional target space. 

Ref. [23] considered various applications where open or closed string background 
gauge fields or gravitational field were turned on. The open string rolling tachyon 
background gives a nice new nontrivial example to test the master formula (I38p . 
The worldsheet action is nonpolynomial, and the master formula involves two-point 
functions in the interacting theory. We choose the local boundary operator to be the 
exponential, O = exp{iuX }, inserted at@ z = e' lT . Its one-point amplitude is fl2T|) . 



which already is a (space)time integral. So we need to show that the integrand Ai(x°) 
can be rewritten as a total derivative as in (13"8"1) . Our convention for the metric of the 
disk is ds 2 = dwdw, so the relation to check is 



dB(x c 



dx° 



(39) 



where 



B(x°) 



Af 



d w 



disk 



w 2 — e 2tT 



2e* 



TBL 



+ 



w — e 



-2lT 



2e 



-(dX {w,w)e^ x0(T) y q 



TBL 



(40) 



where the primes indicate that we have separated the zero mode x°. To show that 
this relation holds, we evaluate the right-hand side. The details of this calculation 
are relegated to Appendix [A] in part because they involve a step that is discussed in 
the next Section [31 The end result is that (139]) holds, so that the one-point function 
is consistent with the general expectation from (]38"1) . 



3 On Coulomb gas relation 

The TBL is related to a statistical mechanical system, the Dyson gas of particles on 
a unit circle [26,27]@ The key property is that the two-dimensional Green's functions 

7 The one-point function is eventually independent of the location. 

8 The analogy has recently been extended to full S-brane (or timelike boundary sine-Gordon 
theory) [28] and to non-BPS half S-brane [29]. 
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can be interpreted as coming from two interacting Coulomb gas particles confined on 
a circle, 



V{e 



log \e %u 



(41) 



where ti, tj are the respective angles. The perturbation expansion in A of ([3]) becomes 
related to the grand canonical ensemble of unit charges on the circle, 



oo jy 
Z 



N=0 



' N 

n 

i=l 



2tt 



-/3H 



(42) 



where the inverse temperature is fixed to (3 = 2, z is the fugacity, and N\ accounts 
for identical particles. The Hamiltonian contains only a potential energy teru0 

h = J2 v(u, tj) = - J2 lo s i e * - ^ i • ( 43 ) 

pairs l<i<j<N 

In this paper we focus only on the canonical ensemble. As discussed in [27], 
correlators in TBL are related to adding additional particles into the ensemble. The 
one-point function ([TO]) requires one additional particle with an arbitrary charge £ at 
an angle r. The Hamiltonian becomes 



y~] log \e lU - i "' i - 



H e = - y . io si e 



l<i<j<N 

and the canonical partition function is 



i<i<j<N 



\e" - e 



iU i 



(44) 



dr 

,V! ./ 27 

dr 

N\ I 2^ 



n 



dU 
2tt 



n 



dtj 
2tt 



n 



I it- it ■ I 2 



n 



(45) 



(46) 



The integrand does not depend on the angle r, hence it can be consistently set to 
zero. We recognize = I%(N) of ffTTj) . 

We can now draw insight from the physical interpretation to better understand 
the integrals and their various extensions. As an example, consider the integral 
corresponding to the canonical ensemble expectation value 



A? 



^cos(r - ti) 

i=l 



1 1 



e lti - e ltj i 



n 



e — e 



i<j 



(47) 

which corresponds to the sum of the projected relative distances of the original charges 
to the additional charge. In part by inspired guesswork we have found a result 



N 



cos(r — ti 



i=l 



(Ah 



3 See the discussion on the physical interpretation in the original paper by F. Dyson [26]. 
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for the integral. We have not constructed a proof for this formula, but have checked 
it for £ = 0,1,2,3,4 and for any N, and a consistency check will be given in Ap- 
pendix [B]0 We can visualize the £ — > oo limit (at finite N) of the result (1481) easily in 
Figure [2j as the additional charge becomes stronger, it forces the unit charges further 
towards the antipodal point of the circle. 



Figure 2: Depicted is the interpretation of (|48p . On the unit circle, embedded in a 
heat bath, there are iV positive unit charges and an additional positive charge £. As 
the charge strength £ increases, the repulsive force acting on the unit charges wins 
over their mutual repulsion, forcing the unit charges closer to each other on the other 
side of the circle. 



4 The n-point boundary amplitude 

The full n-point amplitude ([7]) is very complicated, and so are the integral coefficients 
(Q even at small N, n (> 1). In this Section we will consider an approximation 
or a toy model version of a full calculation. We begin by studying the integrals 
0. We interpret them as Toeplitz determinants. One can then consider a known 
approximation in the large N limit, and try to improve it to be good enough to be 
used in the series expansion (j7]) at every N, while hoping for it to be simple enough 
so that the series can be summed. We use the Coulomb gas analogue and find a 
physically motivated improved asymptotic approximation of (Q. This approximation 
agrees with the previously known asymptotics at leading order in 1/N, but reproduces 
the next-to-leading 1/N corrections to the asymptotics of the integrals better than the 
old result (but not exactly). Even more importantly, it is found to work well for small 
values of N even up to iV = 0, which contribute significantly in the final amplitude 
in the end. However, the approximation is still simple enough to sum the series in 
iV to calculate the integrated n-point amplitude. In the approximation, essentially 

10 After the first version of this work was finished, we were informed by H. Schomerus that he has 
constructed a proof [30] of this formula. We thank him for bringing this to our attention. 
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the "interactions" between the £ a insertions can be neglected. In our end result, the 
n-point amplitude factorizes to a product of n independent one-point amplitudes. We 
also present a simple example that helps to understand and motivate the derivation 
in Appendix iBl 



4.1 Large N asymptotics 

To start with, the t{ integrals in (Q can be done [15, 16] giving 

dr. 



.(AO 



n 

a=l 



2tt 



n 



X<a<b<n 



det T N [f] 



(49) 



where detT/v[/] is the N x N Toeplitz determinant of Fourier coefficients of the 
function 



/n,.,r B (*)=ni e * r- - eft | a6 " > 



(50) 



0=1 



see [15, 16] for more details. 

The determinant is too complicated to allow us to sum the series (J7|). However, 
Toeplitz determinants are known to simplify at large N. In particular, the large A" 
asymptotics of the determinant detTjv[/] is known for fl50l . It reads [31,32] (see 
also [33]) 



detT N [f] = N^iZ JJ \e*--e*\-™fl %% + ^ 

l<a<b<n «=1 G ^ a + ^ 

Moreover, the asymptotic behavior of (l4"9~I) factorizes, 
7«, £ J.V) = detT N [f] H \ e ** - e*»\*>*> 



l + O 



N 



(51) 



o=l 



l<a<b<n 
G(£ a + l) 2 

G(2U + 1) 



l + O 



N 



A G(£ a + l) 2 G(N + 2j a + 1)G(N + 1) 



l + O 



N 



a=l 



1 + °']v 



(52) 



where T^ lr ^ n (N) is the asymptotically r a independent factor of detT N [f] and we 
used (!20l) to write the asymptotics in terms of Barnes G functions (see also (|26|) ). 
Here T^(N) = I^(N) is the one-point function discussed above in Section [2j 

The asymptotic formula (1521) has a nice physical interpretation in terms of the 
classical Coulomb gas on a circle, where ,,..,£„ (A") is the partition function for A" 
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identical unit charges at the inverse temperature j3 = 2, with n additional particles 
having charges £1, . . . , £ n a t fixed angles t%, . . . , r n . Let us assume for a moment that 
all £ a are positive integers^ Then each particle with charge £ a can be thought to be 
a cluster of £ a unit charges. We can then imagine constructing a typical configuration 
of the gas with the n test charges, from a gas of N + un ^ charges, by clustering 

unit charges at distinct locations to form the test charges £ a . For n < Yla=i £a ^ ^, 
the typical separation of unit charges is ~ 1/iV, much less than the typical separation 
between the test charges/charge clusters. Now, we can first interpret the factors 
in (1521) arising from the self-energies of the charge clusters. For a cluster with charge 
£ a , the self-energy is given bjo 

E self = - teg\xi-x 3 \ ~ Yl log iV ~ log N (53) 

l<«<j<£a l<«<j<£a 



giving the contribution 



e 



« ^ ^ (54) 



to the partition function. Second, the factorization of fl52|) can be understood as 
the absence of intercluster interactions at this level of approximation. A heuristic 
argument could be the following. Consider a large number of unit charges on the real 
axis (a piece of the unit circle after magnification) with a typical separation d ~ 2ir/N . 
Choose ( C iV charges at x±, . . . ,x% near the origin (so that x^ ~ d) and perturb 
their locations, Xk — > Xk — 5k, by 5k ~ d symmetrically such that Yl\=i &k — 0, to 
create a cluster of charge £. The change in the electrostatic potential after creating 
the cluster is then 

c 

5k 



1 + 



X Xk 



(55) 



AV(x) = -J^log 

k=l 

as felt at point x outside the cluster, x > Xk- For x ^> d we find that the deformation 
of the potential vanishes rapidly, 

cT 2 



and the contribution to the change in the total energy in the leading order must thus 
come from the interaction of the cluster between the unit charges within the region 
x ~ d ~ 1/N. However, at the distance to the neighboring clusters, the change is 
negligible, so the intercluster interactions are suppressed. 



n This is potentially a dangerous assumption, since eventually we want to set £ a = iu> a where typ- 
ically uj a are real, and naive continuation from integers to complex plane is known to be problematic 
- see, e.g., the discussion in [15]. We will return to this issue in the end of the Section. 

12 In fact, one obtains E se \f ~ £ a (£a — 1) log N/2, but the term — £ a logiV/2 cancels against a change 
in the 1/iV! factors which are discussed below. 
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The analysis suggests a natural way to try to improve the asymptotic formula 
As explained above, in (1521) intercluster interactions are absent. However, the 
1/N corrections due to these interactions can be added almost completely in a very 
simple manner. Naturally, each of the £ a charges must feel the Coulomb force of the 
iV unit charges and the other clusters with £& charges, a ^ b. An easy modification 
of (|52j) to accommodate these is the following. Increase the number of background 
unit charges acting on the £ a cluster from N to N a = N + J2b^a £& m ^ ne asymptotic 
formula (152]) . This indeed replaces the Coulomb force of each charge as a cluster of 
separate unit charges. The total effective number of unit charges in the gas then 
becomes N = N a +^ a = N+J2a=i Ca- We write the improved asymptotic formula for a 
renormalized Toeplitz determinant T^ 1} _^ n (N), which is simply related to T^ 1} „^ n (N) 
of (1521) . The modification is needed since T^ lt ^ n {N) contains the normalization factor 
1/N\ which we want to be replaced by 1/N\: 



i! £ 



.(AO 



r(jv + i; 
i 

r(iv + i) 

" N n 

nn 



-Tt 



,(N) 



(57) 



n 



l<a<b<n 



N 

n 



2tt 



n ^ 

l<i<j<N 



Ati\2 



=1 a=l 



Following the discussion above, we replace the asymptotic formula (1521) by an 
improved formula for ([Sj 



a=l 



where 



4 (Na) 



T(N a + 1) 
r(A>+l) 



F(N+1) 



N a 

n 



dL 



i=l 



2ty 



1 + 



I g^- 7 I ^ 



N„ 



(5? 



JJll-e*!^- 



i=l 



(59) 



n 

l<i<j<N a 

is the properly normalized partition function for a £ a charge in the background of 
N a unit charges. Inverting the relation (1571) . we can rewrite (|58|) as an improved 
approximation for T^ u ...^ n , 



T^ n (N) » T m l[T^(N a ) 



(60) 



a=l 



rp n G(N + £ a + 1)G(N -Ca + i) G(U + l) 2 _ 
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where we introduced the notation T^ rx ^(N) for the improved asymptotics and the 
normalization factor reads 



T 

-L Til 



r(N + r 



n-l 



(61) 



Note that T^ IX c reduces to fl52l) for N — > oo and still has 1/N corrections, but they 
are expected to be essentially smaller than for ( 1521) . In Appendix [B] we discuss the 
simplest nontrivial example (£1, £2) = (2, 2), where the exact results are known [34,35] 
and find that the improved asymptotics AM]) reduces the deviation from the exact 
result by more than an order of magnitude at large N. The new asymptotics continues 
to be a very good approximation to the exact result even for small values of N. 
Moreover, note that setting, e.g., £ n = 1 in (1B~4"j) correctly reproduces T^ rx ^^(N + 1) . 

Finally, we collect our results in a new asymptotic approximation for the Toeplitz 
determinant: 



detT„[/] 



n 

~a<b- 
n 

n 



l<a<b<n 
n 

X 



r(jv + i) 



x 



„_, r(iv + E 6 ^ + i) G(2e + i) 

g(jv + & + E b 6, + i)g(jv - £» + E b 6 + 1) 
c(iv + E & 6 + i) 2 



(62) 



If we substitute this to (1491) . we note that the integrals over r a give 



'6=o,...,£ 



„=o(iV) 



nar a 
~2k 



a=l 



n 



,K»<Kn 



AT 



(k a ) , (63) 



and the result for the integral becomes 



6.W 



raprx 



(iV) 



AA (fc fl ; 



r(iv + E a e a + i) 



X 



X 



n 



r(iv 



r(N + i) 

ea + E fc 6 + i)G(e a + i) 2 



GjN + e a + E b 6 + i)g(jv - Ca + + 1) 
G(iv + E fe 6 + i) 2 



(64) 



The fcf dependence thus completely factorizes into the normalization factor Af '. 

There is, however, a caveat in the above derivation: the result (1641) only makes 
sense when the normalization integral M is convergent. From the definition (163]) we 
see that the integral is singular whenever any of the products k a ■ hi, — > — 1/2, which 
can easily occur for physical momentum values. These singularities are unphysical and 
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they are absent in the original integral of ([9]). What happens is that for k a -kb — > —1/2 
the t integrals become heavily peaked at r a ~ t^. More precisely, the dominant 
contribution to the integral comes from the region where r& — r a ~ 1/JV. In this 
region the large iV limit does not reproduce the r dependence correctly: for positive 
£ a £b the integrand vanishes more rapidly than suggested by the large N limit as 
T a — T b 0, which creates a cutoff for the normalization integral M . 

To avoid this caveat we shall assume that k a ■ kj, > —1/2 which can be satisfied 
together with momentum conservation only if all spatial momenta are small. Even 
when this condition is not met, the result ([Mil may work as a reasonable model for 
the £ and N dependencies of (|4"9"1) . We are planning to study the r dependence of the 
Toeplitz determinant more closely in a forthcoming publication. 

4.2 A model amplitude 

Let us now study what can be said about the integrated amplitude 

A n (2nXe x °) (65) 

oo 
N=0 

based on the asymptotic formula (1641) . Notice that since the coefficients (N) are 

asymptotically equal to a product of one-point functions they also exhibit a powerlike 
behavior for large N [see f)26p ]. This fact strongly suggests that the analysis of 
Subsection 12.21 can be extended to higher point functions, which requires that there is 
an analytic continuation (s) of to complex values of N = — s that 

has a powerlike behavior for s —>■ oo in all sectors of the complex s plane. At least, 
as we shall see below, the continuation exists for the asymptotic formula (|64p (and 
also flo"2"j) ). Also, we calculated I^ lr ^ n (N) for sets of small positive integers £ a and 
for k a ■ kb = in [34], and found that they are, in fact, polynomials of N. See also 
Appendix IB1 where we treat a simple case, (£1,^2) = (2,2), as an example. 

This motivates us to check what is the result if one simply inserts the improved 
asymptotic formula (IMI) to d7j) and to repeat the analysis of Subsection 12.21 For the 
sake of concreteness, we discuss the two-point functional. As explained above, it is 
required that there is such an analytic continuation of I^ V ^(N) of fl64|) to complex 
s = —N that does not blow up exponentially for \s\ — > 00. Remarkably, the simplest 

13 For the two-point function spatial momentum conservation and on-shell conditions actually fix 



A 

sin 



dx° exp 



x 



dx c 



exp 



x 



a=l 
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continuation of ([Mil works: 



JTTfa 



2 



n G{j a + 1) 2 g(-g +ei+6+eo+ i)G(2 - g + , , 

^(2^ + 1) G (_ s + 6 + 6 + i)2 ^ J 

indeed has a powerlike asymptotic behavior for large \s\ as can be verified using the 
formulae flMJ and (|52p above. Further, we need to check that the singularities of 
J^h{s) do not conflict with the contour deformations of Subsection 12. 2[ If £ 1)2 > 
the poles of J^^(s) are located at s = £i +£2 + 1, £i +£2 + 2, . . .. As for the one-point 
amplitude in Subsection l2.2[ they are to the right of s = £i+£2, where we will evaluate 
Jfrfa(s) in the end [see ([68]) below]. As discussed in Subsection 12.21 this means that 
the model two-point function 

00 

AT x (z) = E(-^^( JV ) ( 67 ) 

N=0 

has no singularities for z > and vanishes sufficiently fast for x° — > 00 to make the 
integral over x° in ([65]) convergent. Notice that this is not the case if the "naive" 
asymptotic formula fl52l) is used instead of ([64]) 1^1 In particular, as discussed in 
Appendix [Bj ( [ST]) has the correct asymptotic behavior for z — > +00 for integer (£1, £2) 
only if one uses (IB"!]) . 

The above checks ensure that following the analysis in Subsection 12.21 (see equa- 
tions (151]) . (|35|) . and ([371) ) we can sum and integrate the approximated integrals if^* 
of (161]) . The result is a model two-point amplitude, 



^ 2 W Wfe sinh 7 r(c 1 + c 2 ) + 



a=l 



where we already rotated to imaginary £ a = 2c<j a and omitted a 5-term. Notice the 
similarity to f|2Tj) which stems from the factorized form of the asymptotics (161]) . 

Similarly as the one-point amplitude in (ETJ) , the final result ([68]) is expected to 
include a term oc The delta term arises in the x° integration of A.2{x°) from the 

oscillations in the region x° — > —00: indeed, for imaginary £ = £1 + £2 = + iu 2 the 
integrand e 1 ^ 1 ^ 2 ^ A2(x°) continues to vanish exponentially for x° — > +00, but for 



14 However, also the naive formula leads to a well-defined integral for imaginary £ Q = iuj a which 
we will need in the end. 
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x° — > — oo the function „4 2 approaches a constant, which leads to oscillating behavior. 
The resulting 5 contribution can be isolated as follows. Write 

A 2 = S 0;M2 J dx^+^AT^x ) 

+ / dx°e^ + ^° [Af^(x°) - I^ 2 (N = 0)Mx°)lJ | (69) 

where Ai(x°) |^_ = Ao(x°) = l/(l+27rAe a:0 ). Then the integrand of the first term has 
a simple form and oscillates for x° — > — oo while that of the second one is complicated 
but vanishes exponentially in both directions x° — > ±00. Hence the 5 contribution 
comes solely from the first term which can be integrated exactly, while for the second 
integral is well defined even for imaginary £1 + ^2 and the analytic continuation of 
(168]) from the region of Re(£i + £2) > can be trusted. The 5-term that adds to (1681) 
is seen to b«0 

A* = ^^JZJ^N = 0)5(^1 + uj 2 ) . (70) 

Naively, since we effectively replace iV — >• —i(ui +U2) in the asymptotic expansion 
one would expect the improved asymptotic formula (J55|) to be a good estimate 
for large energies for which 1/N is small. Unfortunately, the correction to AMI) is 
likely to include terms which are ~ uo a /N (or worse) and become (at least) 0(1) at 
N ~ u a . Thus the result ( J68l) only serves as a model of the exact result for any values 
of the energies. However, note that the calculations of Appendix [B] suggest that the 
correction terms are small: the improved approximation is seen to work well also for 
small values of iV and also when continued to negative (but small) values of N. 

It is straightforward to check that the analysis of Subsection 12.21 can be similarly 
applied to the approximate asymptotic formula (1621) for n > 2. The consequent 
generalization of (1651) can be simplified to read 

a. 6 a ^u [<soi (^)-^T (>5>) n G "^^r^ ) ■ < 71 > 



a=l 



As discussed in Subsection 14.11 the k a dependence factorizes into the factor Af (k a 

at least for small k a . We give the singularities and the asymptotic behavior of the 
model amplitude in Appendix O 

Notice that our approximation is not restricted to only positive integer valued £ . 
The original asymptotic formula ( 1321 is valid (see [33]) for Re£ a > —1/2. While our 



15 Naturally, a similar term also appears in the exact amplitude, and is given by eliminating the 
superscript "aprx" of I in ijTD]). 



21 



improved formula (I64p was motivated using integer £ a , it approaches the original one 
(1512"]) at iV — > oo for any set of (complex) £ 's. In Appendix|B]it was demonstrated that 
the improved asymptotics works much better than (I52p for sets of integer £ a . There 
is no reason to believe why it should fail to be an improvement also for imaginary 

Final comments. We conclude with some final thoughts. We have presented a 
method to calculate n-point boundary functions. It would be important to develop 
similar methods for n-point bulk correlators. The main obstacle for a straightforward 
generalization of our calculations is the following. The boundary operators correspond 
to test charges that we constructed from unit charges of the Dyson gas. However, 
the bulk operators cannot similarly be made of the unit charges on the boundary - a 
different trick must be found for the bulk correlation function calculations. Another 
important issue is to develop a clear estimate how good an estimate (190]) is for the 
amplitude. A promising way to test our method would be to work directly in space- 
like boundary Liouville theory, use our method to compute the boundary two-point 
function, and then compare with the exact known result of [17-19]. 
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In this Appendix we check ( 155]) . Following [23], equation (4.7), we need to calculate 



A On KRS relation 



B 





iwX°(r) e -S hdly y 



+ 



(dX°(w,w)e 



•x°(r) -s hdiy y 



(72) 



2e~ iT 
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where the primes of the expectation values indicate that the zero mode x° is left 
unintegrated. We start from (straightforward use of Wick theorem) 



N 



c 



o iuj c X (w,w) 



(73) 



i=i 

|1 - ww\-^ 2 \w - e"\' 2wu>c \{ \e* - e^fJI < 1 ■ " |2 '~' 1 • ■ " ^ 

i<j i 



\w — e 



which is to be integrated over ti = . . . 2ir and summed over N. Notice that 

d 2 



(dX°(w,w)e^l[e x0 ^y 



dwdtu, 



-C 



Ll) c =0 



1 £ 
> ^ + — — 

' iy — e 4 — e l 



n 



\ e itt _ e it 3 \ 2 JJ\ e iT -e iu \ 2 ^ 



(74) 



and similarly for the term containing d in (1721) . Recall that £ = io;. Let us do the w 
integration first. The ^-dependent part reads 



d w 



w 2 — e 2lT 



2e l 



1,, _ 



+ 



+ h.c. 



w — e' H u> — e' 



(75) 



where the h.c. assumes real £. Developing the integrand at w, w = we see that only 
the constant term survives, 



Iw = + f ( eiT ~ ik + e~ iT+lk ) 



7T 



£ + cos(r — £j 



(76) 



Let us then do the integrations. The integral oc £ is the I(-(N) discussed above, and 
for the second term we use (j4"5|) . Putting the results together, 



B 



r2 /■« 



\2V 



IT 



N 



(77) 



o 



Moreover, the x° dependencies of the l/zfi and the e^ x ° exactly cancel, whence after 
derivating with respect to a; [23] 



dB 

dx° 



(7? 



We have thus checked the formula (3.14) of [23] in this special case. The £ 2 in the 
proportionality constant arises from the conformal dimension of the operator e^ x ° , 
which is included in the normalization factor Af of Subsection 12.31 
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B A special case of the n-point amplitude 
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To clarify the involved derivation of the model for the n-point amplitude in Section 
we consider here the simplest nontrivial example, (£1,^2) = (2,2) and k\ ■ k 2 = 
that can be calculated also exactly. We have also checked other cases of sets of small 
integers, and found similar results. 

Let us start with the results for the integral I^iN) of ([9]) which equals h(N, 2)/N\ 
of [34]. Hence we have 



h,2{N) 



35 



(N + 8)\ (N + 7)\ (iV + 6)! (iV + 5)! 
M (JV-1)! (JV-2)! (JV-3)! 



35iV 3 + 467iV 2 + 2046iV + 2940 
5040 



k=l 



,-asymp 
i 2,2 
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(N + 3)(N + 2) 2 (N +1) 
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(N + 2)(N + 1) 
(iV + 4)(iV + 3) 

(N + 5)(N + 4f 
144 



(N + 3) 2 (N + 2) 4 (N + 1)< 



(N + 5i)(N + 4) 2 (iV + 3) 



144 
-1 2 



12 



H(N+k) 



(79) 
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where the first, the second, and third formula are the exact result, the one obtained 
from the asymptotic formula (1521) . and the improved asymptotic formula flM|) . respec- 
tively. For iV — > 00 we find 



y-asymp 
i 2,2 



(N) 



432 142384 



h,2{N) 
(N) 



Taprx 
i 2,2 



h,2( N ) 
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35A^ 
12 
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1225iV 2 
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1 
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Some of the values of the integrals are tabulated in Table [U The improved formula 
works much better, in particular for low values of N. For higher values of £ a the 
improvement is even more drastic, basically since the difference between the effective 
number of unit charges N = N + ^2^=i £o> which is used in the improved asymptotics, 
and the actual number of unit charges N increases. 
The two-point function 



A 2 (x°) 



1*1=6=2 



N=0 



16 Notice that the condition k± ■ = eliminates all dependence on spatial momentum. It actually 
conflicts with momentum conservation. 
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Table 1: Exact and approximated values of the integral 22,2 (2V)- The first tabulated 
row is the exact result of the integral, while the two others are given by the asymptotic 
formulae (152|) . (1641) . written explicitly in ([79]). 



N 





1 


2 


3 


5 


10 


100 


h,2(N) 

2 2 T P (iV) 


70 
1 

200/3 


784 
36 
750 


4590 
400 
4410 


18968 
2500 
54880/3 


175320 

38416 

170100 


7514650 
2944656 
7357350 


91680976745020 
81349594398801 
91384995374400 



where z = 2ir\e x , can be calculated explicitly for all the results ( 1791) . In particular, 
for large x° we have 
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(2ttA) 
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The improved asymptotic formula produces also the large x° asymptotics nicely: A. 
is correct up to the proportionality constant for x° — > oo. 

The analytic continuation of I 2j2 (N), J2, 2 (s), is found by letting iV — > — s in (179|) . 
The result for the integrated amplitude is then obtained by applying (1681) to the three 
different cases of (1791 . which gives 
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(83) 



where the first, the second, and the third numbers are the exact result, the result 
for the naive asymptotic formula (152|) . and the result for the improved formula (l64j) . 
respectively. Both the asymptotic approximations give an incorrect result by an 
infinite factor. However, the numerical factor —1/70 of the exact result is extremely 
small when compared, e.g., to the series coefficients of Tabled! whence the zero result 
obtained by the improved asymptotic formula should be, in fact, considered as a good 
approximation. It would be interesting to be able to compare our model amplitude to 
the exact result for more physical, noninteger values of £ a where no accidental zeroes 
or infinities are expected to occur. 

We end this Appendix by an encouraging observation in a more general setup. It 
is straightforward to check that, in fact, for any sets of integer (£i, . . . £ n ) the expected 
exact asymptotic behavior [34] 
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is reproduced by *4^ prx similarly as for £ 2 ) = (2, 2) in (|82|) . We denote the analytic 
continuation of I^ rx ^(N) of to complex s = — iV by J&™£ { s ) i n analogue to 
( TT8l) . (1661) above. Using the behavior of Barnes G near its zeroes from Appendix O, a 
lengthy calculation shows that the first pole of ^|f rX f n ( s )/ sin7rs on the positive real 
axis occurs at s — J2 a £a + ma x a {£ a } in the special case of integer £ a . Hence, for the 
n-point function and integer £ a , (15Uj) indeed becomes 



AT x (x°) ~ exp 



x°^2^ a - x° max{£ a } 



55) 



In other words, the corresponding poles of J r |f >rX £ n and the exact analytic continuation 
J&,..4n coincide. Note that these poles lie at positive s, i.e., negative N, while J^^ n 
results from an asymptotic formula (jMj) for large positive N. This observation gives 
more confidence to the model amplitude of (1681 . (ITT]) which was derived using (1641) . 

C Singularities and asymptotics of the model am- 
plitude 

Barnes G(z) is an entire function and has its zeroes on the negative real axis, 

G(z + 1) = {-l) n{n ~ 1)/2 G{n + l)(z + n) n [1 + O {z + n)] , (86) 

for n = 1, 2, Hence all the special points (zeroes or singularities) of the one-point 

amplitude A\{oS) of ff2~Tj) lie on the imaginary axis, 

A = ^ - {iuj-n) [l + 0(iu-n)\ ; 
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for any n = 0, 1, 2, In particular, poles are found at u = (where Ai ~ l/iuf) 

and at u; = i/2, 3i/2, hi/2, . . .. 

The singularities of the model n-point amplitude 

\ a / a=l v ' 

arise similarly from the poles of V (i uo a ) and from the zeroes of each G(2iu a + 1). 
As above, the hat denotes that we dropped the k a dependent terms. At the possible 
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singularities Ea 1 ^ — —m, iujb — —m, iuJb — m+ 1, and iuJb — —m — 1/2, we find 
(-27rA) m -A- G{iu a + 1) 3 G(2 - iu a ) 
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respectively, for each m = 0,1,2,... and b = 1,2, ... ,n, and assuming that the 
singularities are distinct (which requires n > 1). We thus find poles only at negative 
integer values of i E a w a and at negative half-integer values of each iuj^ whereas the 
amplitude vanishes for almost all integer values of each iui,. 

Since the model amplitude has a factorized form, its asymptotics can be immedi- 
ately written down based on the one-point result (|25|) . It reads 
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a=l 



Note that the prefactor vanishes exponentially for any uj a — » oo and thus does not 
affect the leading asymptotic behavior ~ expfcu 2 ], which is similar to that of the 
one-point amplitude in 
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